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ABSTRACT 

Using a simple analytic approach we address the question of whether radiative cooling, nongravita- 
tional heating and cooling plus heating models can simultaneously explain the observed global X-ray 
properties (entropy and X-ray luminosity distributions) of groups and clusters and the residual soft X- 
ray background (XRB) after discrete sources are removed. Within the framework of typical cold dark 
matter structure formation characterized by an amplitude of matter power spectrum cts = 0.9, it is 
argued that while radiative cooling alone is able to marginally reproduce the entropy floor detected in 
the central regions of groups and clusters, it is insufficient to account for the steepening of the X-ray 
luminosity - temperature relation for groups and the unresolved soft XRB. A phenomenological preheat- 
ing model, in which either an extra specific energy budget or an entropy floor is added to the hot gas in 
groups and clusters, fails in the recovery of at least one of the X-ray observed features. Finally, the soft 
XRB predicted by our combined model of cooling plus heating exceeds the observational upper limits 
by a factor of ~ 2, if the model is required to reproduce the observed entropy and X-ray luminosity - 
temperature relationships of groups and clusters. Inclusion of the cosmic variation of metallicity and the 
self-absorption of the cooled gas as a result of radiative cooling in groups and clusters, or exclusion of 
the contribution of nearby, massive clusters to the XRB does not significantly alter the situation. If the 
discrepancy is not a result of the oversimplification of our analytic models, this implies that either our 
current understanding of the physical processes of the hot gas is still incomplete, or the normalization 
of the present power spectrum has been systematically overestimated. For the latter, both the X-ray 
properties of groups and clusters and the XRB predicted by preheating model and cooling plus heating 
model can be reconciled with the X-ray observations if a lower value of the normalization parameter 
(Ts ~ 0.7 is assumed. 

Subject headings: cosmology: theory — diffuse radiation — galaxies: clusters: general — intergalactic 
medium — X-ray: general 



1. INTRODUCTION 

A substantial fraction of the baryons in the local uni- 
verse exists in the form of diffuse warm-hot intergalac- 
tic medium (IGM) with temperatures of T ~ 10^-10^ K 
as a result of gravitationally driven shocks and adiabatic 
compression as they fall onto large-scale structures and 
collapsed dark halos (Gen & Ostriker 1999). In the for- 
mer case, the warm IGM may escape direct detection, but 
measurements of the strength and power spectra of the 
cosmic soft X-ray background (XRB) (e.g. Soltan, Frey- 
berg & Triimper 2001) and the Sunyaev-Zel'dovich effect 
(e.g. Bond et al. 2002) can set stringent constraints on its 
amount and distribution. For the latter, the hot IGM in 
virialized dark halos such as groups and clusters manifests 
itself as strong X-ray emission sources through thermal 
bremsstrahlung, which provides a powerful tool for study- 
ing not only the physical processes of the hot IGM but 
also the formation and evolution of groups and clusters. 

While gravity and thermal pressure play potentially im- 
portant roles in the overall distribution and evolution of 
the hot IGM in groups and clusters, the IGM also suffers 
from the influence of nongravitational effects such as ra- 
diative cooling, heating by supernovae and/or AGNs, non- 
thermal pressure, etc. In low-mass groups and the central 
regions of clusters, nongravitational effect may even be- 



come dominant for the evolution of IGM. Indeed, there 
is growing observational evidence for the presence of non- 
gravitational effects in groups and clusters. Among many 
arguments, the most convincing observational facts are 
the significant departure of the observed X-ray luminos- 
ity (ix) - temperature (T) relation of groups and clusters 
(Lx (X T^~^) from the prediction (Lx cx: T^) of self-similar 
model (e.g. Edge & Stewart 1991; David et al. 1993; Wu, 
Xue & Fang 1999; Helsdon & Ponman 2000; Xue & Wu 
2000 and references therein) and the entropy [S] excess in 
the central cores of groups and clusters (Ponman, Cannon 
& Navarro 1999; Lloyd-Davies, Ponman & Cannon 2000). 
Two prevailing scenarios suggested thus far, preheating 
and radiative cooling of the hot IGM, both of which tend 
to suppress the X-ray emission of the hot IGM heated by 
purely gravitational shocks and compression, have been 
shown to be indistinguishable in the explanation of the 
observed X-ray properties of groups and clusters (Voit & 
Bryan 2001b; Voit et al. 2002; Borgani et al. 2002 and 
references therein), although each model still has its own 
problem. For instance, preheating model suffers from the 
so-called energy crisis that an unreasonably high efficiency 
of energy injection into the IGM from supernovae must be 
required in order to bring the IGM to the energy level seen 
in the Lx-T and S-T relations of groups and clusters (Wu, 
Fabian & Nulsen 1998, 2000; Tozzi 2001; etc.), although 
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energy supply by AGNs may help to reduce the discrep- 
ancy (Valageas & Silk 1999). Radiative cooling model suf- 
fers from the overcooling problem, i.e., the material that 
has cooled out of the hot IGM greatly exceeds the obser- 
vational limits (Balogh et al. 2001). Moreover, cooling 
efficiency is also a major concern in the explanation of the 
X-ray observed properties of groups and clusters (Bower 
et al. 2001). Nonetheless, both (pre)heating and cooling 
are two natural processes during the formation of galax- 
ies, which is justified by the measurements of the IGM 
enrichment at high redshifts and the stellar mass fraction 
(~ 10%) in the local universe. Indeed, it has been realized 
recently that only can a combination of preheating by su- 
pernovae and radiative cooling of the IGM reproduce the 
observed X-ray properties of groups and clusters (Voit et 
al. 2002; Borgani et al. 2002). 

An independent and sensitive probe of the physical 
properties of the IGM is through study of the soft XRB. 
It has been shown that the XRB produced by the gravi- 
tationally heated and bound IGM in groups and clusters 
within the standard framework of hierarchical formation 
of structures vastly exceeds the upper limits set by current 
X-ray observations, lending further support to nongravi- 
tational heating scenario (Pen 1999; Wu, Fabian & Nulsen 
2001; Bryan & Voit 2001). Yet, it is also noted that some 
hydrodynamical simulations have yielded a soft XRB com- 
patible with current limits without inclusion of nongrav- 
itational heating (Croft et al. 2001; Dave et al. 2001; 
Phillips, Ostriker & Gen 2001). This partially reflects the 
difficulty in handling the complex processes of IGM evolu- 
tion in groups and clusters. Essentially, an excess energy 
of about 1 keV/particle should be injected into the IGM 
in order to reduce the predicted soft XRB to a reasonable 
level. On the other hand, it is unlikely that radiative cool- 
ing alone is sufficient to eliminate the discrepancy (Wu et 
al. 2001). 

The question we would like to address in this paper 
is: Can radiative cooling and nongravitational heating ex- 
plain simultaneously the observed global X-ray properties 
of groups and clusters and the unresolved cosmic XRB ? 
Previous studies have reached the following conclusions: 
If one leaves the energy source problem aside, preheat- 
ing model is in principle able to reproduce the observed 
Lx-T relation and excess entropy in the central regions of 
groups and clusters (e.g. Cavaliere, Menci & Tozzi 1997, 
1998; Balogh et al. 1999; Tozzi & Norman 2001; Babul 
et al. 2002; etc.) and predict a consistent XRB with that 
observed (Pen 1999; Wu et al. 2001; Bryan k Voit 2001). 
As for radiative cooling, it has successfully accounted for 
the entropy floor seen in the central cores of groups and 
clusters (Voit & Bryan 2001b). By properly truncating the 
outer radii of X-ray surface brightness for groups and clus- 
ters in terms of current X-ray flux limits, one is also able 
to reproduce the observed steepening of the Lx-T relation 
of groups and clusters (Wu & Xue 2002a; Voit et al. 2002). 
In particular, the Lx-T relation found by hydrodynamical 
simulations in terms of cooling shows a good agreement 
with observations (e.g. Muanwong et al. 2001, 2002; Bor- 
gani et al. 2002). Consequently, it is crucial to reexamine 
whether or not cooling model is also sufficient enough to 
suppress the contribution of the IGM to the soft XRB. If 
the negative result of Wu et al. (2001) is confirmed, it 



deserves to explore the combined effect of radiative cool- 
ing plus heating by supernovae on the X-ray properties 
of groups and clusters and the soft XRB. Failure of all 
these efforts may indicate that our current picture of the 
physical processes of the hot IGM in groups and clusters 
is incomplete at some level. Throughout this paper we 
assume a flat cosmological model (ACDM) of flu = 0.35, 
Oa = 0.65 and h = 0.65. 

2. DARK HALOS 

Dark matter distribution in a virialized halo is assumed 
to follow the universal density profile suggested by numer- 
ical simulations (Navarro, Prenk & White 1997; NFW) 

l^chPcrit 
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(r/r,)(l-f r/r,)2' 

where Sch and arc the characteristic density and length 
of the halo, respectively, and pcrit is the critical density 
of the universe. In order to fix the two free parameters, 
Sch and Vg, we first specify the concentration parameter 
c = J'vir / for a given halo of mass M through the empiri- 
cal fitting formula found by numerical simulations (Bullock 
et al. 2001) 
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Next, we define the virial mass M such that within the 
virial radius rvir the mean mass density of the dark halo 
is Ac times the critical density of the universe: 



M = -Trr^AcPcrit, 



(3) 



where for a flat, ACDM cosmological model, Ac = IStt^ -|- 
S2[nM{z) - 1] - 39[17m(^) - 1]', ^m{z) = ^Mil + zf/E^ 
and E'^ = flui'^- + z)^ + Finally, we determine the 
virial temperature using cosmic virial theorem (Bryan and 
Norman 1998): 
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In which we have taken the normalization factor to be 
/t = 1. We have also tested a lower value of /t = 0.8, 
and found that our results remain almost unchanged. 

We use the modified PS mass function by Sheth & Tor- 
men (1999) to describe the abundance and evolution of 
virialized dark halos that grow from random-phase Gaus- 
sian initial fluctuations: 
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where A — 0.3222, a = 0.707, p = 0.3, p is the mean cos- 
mic density, 6c is the linear over-density of perturbations 
that collapsed and virialized at redshift 2;, cr is the linear 
theory variance of the mass density fluctuation in sphere 
of mass M = Att/jR^ /3. We parameterize the power spec- 
trum of fluctuation P{k) oc k"T'^{k) and take the fit given 
by Bardeen et al. (1986) for the transfer function of adia- 
batic CDM model T{k). The primordial power spectrum 
is assumed to be the Harrison-Zel'dovich case n = 1. The 
mass variance for a given P{k) is simply 
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where W{x) = 3(sin.x — xcosx)/x^ is the Fourier repre- 
sentation of the window function. The amphtude in the 
power spectrum is determined using the rms fluctuation on 
an 8 Mpc scale, as- We adopt a normalization param- 
eter of (T8 = 0.9 suggested by weak gravitational Icnsing 
measurements (see Refregier, Rhodes & Grith 2002 for a 
recent summary), and then demonstrate the influence of ag 
on the evaluation of XRB using a lower value of (Tg = 0.7. 

3. IGM 
3.1. Global quantities 
Following conventional definition, we use 

~ (7) 
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to measure the entropy of IGM, in which T and rie are 
the electron temperature and number density, respec- 
tively. The total X-ray luminosity in terms of thermal 
bremsstrahlung is 



= y n,nHA{T)dV, 



(8) 



where nn is the number density of hydrogen, and A(r) 
is the cooling function which is calculated from the 
Raymond-Smith (1977) code. Whenever theoretical pre- 
dictions are compared with observations, we work with the 
emission-weighted temperature 

^ jT{r)A{T)n,nHdV 

jA{T)n,nudV ' ^ ' 

Furthermore, we assume that the IGM with and with- 
out cooling/preheating is always in hydrostatic equi- 
librium with the underlying gravitational potential of 
groups/clusters dominated by dark matter: 

1 d[n,{r)kT{r)] _ GMuuir) 

Iimpnc{r) dr ' 

where /i « 0.593 is the mean molecular weight. Finally, 
we evaluate the total XRB intensity at frequency u by in- 
tegrating the X-ray emission of all halos over mass range 
and redshift space 
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where /ip is the Planck constant, Dl is the luminosity dis- 
tance, and dN/ dVdM is the mass function given by equa- 
tion (5). 

3.2. No cooling and no heating 

In the absence of radiative cooling and nongravitational 
heating, we assume that IGM traces dark matter 

Pgas(?') = /bPDM(r). (12) 

We solve the equation of hydrostatic equilibrium under the 
boundary restriction T(oo) = 0, which yields 



kT{r) = kT*— {l + 



vji: 



[l + x] ln(l -hx) - X 
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dx. 



(13) 

where kT* = AttG p^m^ScYiPc-citfl ■ This allows us to calcu- 
late straightforwardly the X-ray luminosity, entropy and 
total XRB. Our numerical computations show that these 



quantities are rather insensitive to the radial variation of 
temperature described by equation (13). Consequently, 
we may adopt a constant temperature instead of equation 
(13) in the evaluation of the global X-ray properties of 
groups and clusters and the XRB. It is well known that 
the results predicted by this so-called 'self-similar' model 
(hereafter Model I) depart remarkably from X-ray obser- 
vations. Here we use this model as a reference point only. 

3.3. Radiative cooling 

Conservation of energy ensures that the energy loss due 
to bremsstrahlung emission is balanced by the decrease in 
the specific energy of IGM in groups/clusters: 



-n^kT = n^nnA{T%, 



(14) 



where rit is the total number density of the IGM. Setting 
the cooling time tc to equal the age of groups/clusters, 
or simply the age of the universe at the redshift of the 

groups/clusters determines the cooling radius TcooI and 
the total mass Mcooi of the cooled material if the IGM 
is assumed to follow the dark matter distribution before 

cooling. Following Voit & Bryan (2001b) and Wu & Xue 
(2002a), we can find the distribution of the remaining IGM 
after cooling by solving the equation of hydrostatic equi- 
librium under the conservation of total baryonic mass 



Mcool(rvir) + Mgas(rvir) = fhM, 

and the conservation of entropy 

T{r) T°[f) 



(15) 



(16) 
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where r is related to f through 

Mga.(r)=MO,,(f)-Meool. (17) 

Here we use superscript ^ to denote the quantities be- 
fore cooling. The above two equations demonstrate that 
the IGM originally distributed between TcooI and f before 
cooling is transported to a smaller region of radius r af- 
ter cooling. We consider an evolving metallicity model of 
Z = 0.3ZQ(t/to) (Model II) and a constant metallicity 
model of Z = 0.3.^0 (Model HI), where to is the present 
age of the universe. 

3.4. Preheating 

A phenomenological treatment of (pre) heating is to raise 
the specific energy or entropy of the IGM in groups and 
clusters to a certain level regardless of whatever the en- 
ergy sources would be. Usually, an energy budget of 0.1 

3 keV/particle or an entropy of 50 400 keV cm^, de- 
pending on the epoch and environment of (pre)heating, is 
needed in order to reproduce the observed X-ray proper- 
ties of groups and clusters. Here we consider two simple 
approaches to demonstrating the effect of (pre) heating on 
the X-ray properties of groups and clusters. 

Model IV: We begin with the IGM distribution pre- 
dicted by self-similar model (Model I), and then simply 
raise the specific energy of each particle in groups and 
clusters by a constant amount of kAT: 

kTir)=kAT + kT,^i{ir). (18) 

The corresponding distribution of electron number density 
ne(r) can be obtained by solving the equation of hydro- 
static equilibrium. A critical point in such an exercise 
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is the boundary condition. To a first approximation we 
may assume that the gas density at virial radius r^ir is 
universal, nel^^vir) = (/b/MeTOp)PDM(''vir)- Nevertheless, 
one should keep in mind that this restriction may fail for 
low-mass systems. Another commonly adopted, inaccu- 
rate method is to take the total gas mass fraction within 
Tvir to be universal Mgas(r"vir) = fbM. Actually, one still 
has no better choice of boundary condition in the deter- 
mination of the gas distribution from the equation of hy- 
drostatic equilibrium. 

Model V: Instead of raising the specific energy of the 
IGM particles, we add a constant entropy floor AS to the 
entropy profile given by self-similar model <S'seif(r) 

S{r) = AS + SMr)- (19) 

We then solve the equation of hydrostatic equilibrium to 
get the electron pressure distribution Pe = n^kT 



Pe{r) = 



1 5/2 



dr 



(20) 



5-3/5(^)^2 

where we have adopted the boundary condition Pe{oo) = 
0. Note that this restriction should not be taken too liter- 
ally because the equation of hydrostatic equilibrium may 
break down beyond rvir- 

3.5. Radiative cooling plus heating 

We first work with the IGM temperature and entropy 
distributions, Tcooi('') and 5cooi('')) predicted by cooling 
(Model II) , and then raise cither the specific energy of each 
particle by a constant amount of kAT (Model IV) or the 
entropy by a constant fioor of AS' (Model V) so that the 
new temperature (Model VI) and entropy profiles (Model 
VII) become 

T(r) = AT + Teooi(r), (21) 

and 



S(r) = AS + Scooi(r), 



(22) 



respectively. In a way similar to Models IV and V, we 
can find the new density profile of the IGM for our com- 
bined model of radiative cooling plus heating. Note that 
our heating model differs from the conventional preheating 
model. For the latter, the IGM is heated before cooling, 
while in our model we address the issue of how large an 
extra energy budget should still be added to the IGM after 
the cooled gas is removed from groups and clusters. 

The parameters of the seven IGM models are summa- 
rized in Table 1. 

4. RESULTS 

Our strategy is as follows: For each model listed in Ta- 
ble 1, we first compute the expected entropy and X-ray 
luminosity distributions of groups and clusters, and then 
compare with X-ray measurements, which constitutes a 
fundamental test for each of the proposed models. For 
the heating models (Model IV and Model V) and cooling 
plus heating models (Model VI and Model VII), this allows 
us to work out the most probable energy budget kAT or 
AS with which the observed entropy and X-ray luminosity 
distributions of groups and clusters can be simultaneously 
accounted for. We then compute the XRB predicted by 
these models, and examine whether their predictions are 
compatible with current observational limits. 



4.1. Entropy distribution 

The entropy at O.lrvir against the emission weighted 
temperature T predicted by all the models is shown in 
Figure 1, together with the updated measurements of 
S(O.lrvir) by Ponman et al. (1999), Lloyd-Davies et al. 
(2000) and Xu, Jin & Wu (2001). Essentially, all mod- 
els except self-similar model can roughly reproduce the 
overall S(0.1rvir)-r distributions. In the preheating model 
(Fig. lb and Fig.lc) and the combined model (Fig. Id) we 
have tested a set of energy input values kAT and AS', in 
attempt to estimate the most probable parameters which 
give the closest fits of both the S-T distribution and the 
Lyi-T relation (see below). In several cases we have also 
shown the curves for two sets of kAT and AS around the 
best-fit values in order to demonstrate how sensitively the 
results depend on the input energy/entropy. While there 
are some differences in the predicted S-T relation among 
various IGM models and for different input kAT or AS 
values especially below temperature T ~ 1 keV, a decisive 
conclusion regarding which model provides the most satis- 
factory explanation of the entropy data cannot be drawn, 
because of sparse data points and associated large uncer- 
tainties. It is emphasized that the entropy profile expected 
from radiative cooling alone (Fig. la) is consistent with ob- 
served data (Voit & Bryan 2001b), although with an ad- 
ditional energy supply by preheating model the predicted 
entropy distribution seems to yield a better fit to the data 
points for low-temperature systems of T < 1 keV (Fig. Id). 



4.2. Lx-T relation 

We now compare the bolometric X-ray luminosity - tem- 
perature relations of groups and clusters predicted by dif- 
ferent IGM models with X-ray observations. For the latter 
we use the catalog of X-ray groups and clusters compiled 
by Wu et al. (1999) and Xue & Wu (2000). The updated 
sample contains 57 groups and 192 clusters whose X-ray 
temperature and luminosity are both available. We have 
converted the X-ray luminosity in the Einstein-de Sitter 
universe into the one in the ACDM cosmological model. 

The Lx-T relations predicted by cooling models II and 
III arc shown in Figure 2. While there is a good agree- 
ment between the predicted and observed Lx-T relations 
on cluster scale T > 4 keV, cooling models become to be 
insufficient to recover the observed data at temperature 
below ~ 4 keV. One possible reason for this discrepancy 
may arise purely from selection effect: The theoretically 
predicted X-ray luminosity accounts for all emission in- 
side virial radii of groups, while current X-ray observa- 
tions have not corrected for lost fiux falling out the de- 
tection aperture for majority of groups with /? parameters 
less than 1/2 (e.g. Ponman et al. 1996; Helsdon & Pon- 
man 2000). For this reason, following Wu & Xue (2002a) 
and Voit et al. (2002) we re-calculate the X-ray luminos- 
ity for Model II by excluding the contributions of the IGM 
in the outer regions of groups and clusters set by a given 
X-ray surface brightness limit Sumit- The expected X-ray 
luminosity of clusters is almost unaffected by this trun- 
cation. However, the X-ray luminosity of groups suffers 
seriously from this selection effect, and the resulting Lx 
drops remarkably at low-temperature end. This may par- 
tially eliminate our concern about the efficiency of cooling 
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Table 1 

Summary of the parameters and results for the IGM models. 



Model 


cooling 


heating 


metallicity (Zq) 


S-T 


Lx-T 


XRB (as = 0.9) 


XRB (ct8 = 0.7) 


I 


no 


no 


0.3{t/to) 


fail 


fail 


fail 


fail 


II 


yes 


no 


0.3{t/to) 


ok 


fail/ok* 


fail 


fail 


III 


yes 


no 


0.3 


ok 


fail/ok* 


fail 


fail 


IV 


no 


yes(AT) 


0.3{t/to) 


ok 


ok 


fail 


ok 


V 


no 


yes(A5) 


0.3{t/to) 


ok 


ok 


fail 


ok 


VI 


yes 


yes(Ar) 


0.3{t/to) 


ok 


ok 


fail 


ok 


VII 


yes 


yes(A5) 


0.3(Vio) 


ok 


ok 


fail 


ok 



*Acceptable if a truncated radius is properly introduced in terms of X-ray surface brightness limit. 



in the explanation of L-^-T relation. Moreover, the L-^-T 
relations of groups and clusters with and without the in- 
clusion of the cosmic variation of metallicity (II and III) 
show little evolution at least out to 2; = 1. 

Figure 3 illustrates the Ly^-T relations predicted by pre- 
heating model (IV and V). Adding an extra energy of 
fcAT = 0.25 kcV to each particle in groups and clus- 
ters roughly reproduces the observed Lx-T relation. Note, 
however, that there is a small but significant vertical shift 
between our predicted and observed data in the tempera- 
ture range 2 < T < 8 keV (Fig. 3a). Wc have then tried 
a larger value of kAT = 0.3 keV, which indeed reduces 
somewhat the difference. A further increase of kAT leads 
the entropy curve in Fig. lb to rise too high to be con- 
sistent with observations. When kAT exceeds the virial 
temperature of a given group, wc assume that the IGM 
cannot be trapped in the system. This is responsible for 
the cutoff of our Lx-T relation at low temperature end. 
For Model V, in which we have raised the entropy by a 
constant floor of AS = 120 keV cm^, the predicted Lx- 
T relation of clusters matches perfectly the observed one. 
For low-temperature groups, the agreement becomes only 
marginal. Increasing slightly the value of AS improves 
the match. But, it meanwhile raises the entropy distri- 
bution (Fig.lc). We have tried a value of AS* = 200 keV 
cm^, and found that the entropy at low-temperature be- 
comes too large to be reconciled with the observed data. 
In a word, it turns out that Model V is marginally ac- 
ceptable. Alternatively, our numerical computations show 
that there is no apparent cosmic evolution of the predicted 
Lx-T relations for both Model IV and Model V, regardless 
of whether or not the metallicity varies with cosmic time 
in terms of 0.3^0(^/^0). 

Because either cooling or preheating alone can repro- 
duce or marginally reproduce the observed Lx-T and 
S'(0.1rvir)-7' relations, it is naturally expected that a com- 
bination of the two mechanisms should be successful in 
the explanation of the two relationships. Indeed, within 
the framework of radiative cooling, an additional energy 
input of kAT = 0.1 keV in Model IV and an entropy in- 
put of AS = 50 keV cm^ in Model IIV both result in the 
X-ray luminosity and entropy distributions which agree 
nicely with X-ray observations (sec Figure 4). We have 
also considered the situation of a constant metallicity of 
0.3^0 and at high redshifts out to 2; = 1, and found that 
the corresponding modifications are only minor. 



4.3. XRB 

A considerably large fraction of soft and hard XRB has 
been resolved into discrete sources (e.g. McHardy et al. 
1998; Hasinger et al. 1998, 2001; Mushotzky et al. 2000; 
Giacconi et al. 2001, 2002; Tozzi et al. 2001; Horn- 
schcmcier et al. 2001; Rosati et al. 2002; Bauer et al. 
2002; etc.). The maximum admitted ranges of the unre- 
solved flux at different energy bands have been summa- 
rized in Wu & Xue (2001). Here we only add a new upper 
limit in the 2-8 keV band from the 1 Ms Chandra ob- 
servation of the Chandra Deep Field North (Cowie et al. 
2002): 0.5 x 10" ergs s"^ cm'^ deg-^. It should be kept 
in mind that the diffuse X-ray emission from some of the 
nearby, bright galaxies, groups and clusters has also been 
included in the current resolved soft XRB, although the 
fractions of these diffuse X-ray sources in the resolved and 
unresolved XRB are still uncertain. In this regard, the 
residual soft XRB after the removal of the discrete sources 
may not be taken to be a very stringent upper limit on the 
contributions of groups and clusters. 

The expected XRB spectra from different IGM models 
in terms of Equation (11) are shown and compared with 
the observational upper limits in Figure 5. Meanwhile, we 
also demonstrate the XRB by summing up the contribu- 
tions of groups and clusters described by their X-ray lumi- 
nosity functions (Wu & Xue 2001). The good agreement 
between the XRB produced by the 'known' population of 
groups and clusters and the current limits indicates that 
the diffuse IGM confined in groups and clusters is prob- 
ably the major source of the unresolved soft XRB. This 
may help to eliminate the above concern that some of the 
bright groups and clusters have been resolved and thus re- 
moved from the residual XRB, and the unresolved XRB 
may not constitute a robust constraint on the diffuse IGM 
of groups and clusters. Actually, nearby bright and mas- 
sive clusters (e.g. z < 0.2 and M > 5 x 1O"M0) only 
make a minor contribution to the total soft XRB (Wu & 
Xue 2001; see also Figure 5). 

It turns out from Figure 5 that the soft XRB predicted 
by the self-similar model (Model I) vastly exceeds the ob- 
servational hmits (Pen 1999; Wu et al. 2001). The huge 
difference of up to two orders of magnitude at E ^ OA keV 
implies that the IGM should have a much shallower distri- 
bution than dark matter especially in low-mass systems. 
Preheating was thus advocated as a potentially important 



6 



mechanism to break the similarity between dark matter 
and IGM. 

We first begin with the radiative cooling scenario. The 
predicted XRB remains almost unchanged if a constant 
mctallicity of 0.3^0 (Model III) is replaced by a time- 
varying quantity of 0.3ZQ{t/to) (Model II). Actually, our 
numerical computations show that this conclusion applies 
to all the models in Table 1. In the hard energy band 
E > 2 keV, the XRB produced by groups and clusters 
within the framework of cooling is well below the observa- 
tional limits. By contrast, in the soft energy band E Ki 0.1- 
2 keV the cooling results are about 2-4 times larger than 
the current upper limits placed on contribution from dif- 
fuse IGM to the XRB. This is consistent with previous 
findings by Wu et al. (2001) based on N-body simulations 
of halo merger trees coupled with semi- analytic models. 
We have also studied the contributions of different halos 
to the total XRB, and foimd that most of the soft XRB is 
produced by groups of mass M ~ lO^^^ Mq. Note that very 
low-mass halos make almost no contribution to the XRB 
because they contain very little hot IGM due to their too 
short cooling time (see Fig. 5a). It appears that although 
cooling provides a more or less reasonable explanation of 
the X-ray luminosity and entropy distributions, it is insuf- 
ficient (by a factor of 2-4) to account for the unresolved 
soft XRB. 

We now turn to preheating model. For Model IV, in 
which the specific energy of each particle in groups and 
clusters is raised by a constant amount of kAT = 0.25 keV, 
we have successfully recovered the observed S'(0.1rvir)-T' 
relation (see Fig. lb) and roughly reproduced the Lx-T re- 
lation (see Fig. 3). Our expected XRB from this model is 
shown in Fig. 5b. Unfortunately, we reach a result very 
similar to the cooling prediction discussed above: Model 
IV fails to efficiently suppress the X-ray emission of groups 
and clusters to a level below the unresolved XRB, and the 
difference in the soft energy band E — 0.1-2 keV is ap- 
proximately by a factor of 2-4. 

For another preheating model V, in which we have raised 
the entropy of IGM by a constant floor AS instead of 
kAT for Model IV, we present the XRB in Figure 5c for 
two choices of AS, 120 and 200 keV cm^, respectively. 
The former nicely reproduces the entropy distribution but 
only marginally explains the Lx-T relation, while the lat- 
ter turns to be successful in the recovery of the Lx-T re- 
lation but results in an overestimate of the entropy dis- 
tribution. Now, in terms of their predicted XRB spectra 
alone, the model with AS = 200 keV cm^ becomes accept- 
able, while the entropy level of AS' = 120 keV cm^ is not 
sufficiently large to reduce the XRB to the observational 
limits. This is roughly consistent with simulation results 
that an entropy floor of 100 — 200 keV cm^ is needed to 
suppress the expected unresolved X-ray background below 
the observational limits (Voit & Bryan 2001a). 

Finally, we come to the cooling plus heating model. 
Adding an energy budget of kAT = 0.1 keV to each parti- 
cle in addition to radiative cooling (Model VI) has allowed 
us to reproduce remarkably well the entropy and X-ray lu- 
minosity distributions of groups and clusters, as are shown 
in Fig. Id and Fig. 4. The XRB predicted from this model is 
consistent with the observational constraints above T ~ 1 
keV (see Fig.5d). Nonetheless, the theoretical prediction 



becomes to be larger than the observational limits by a fac- 
tor of 2-3 in the energy range from 0.1 to 1 keV. Because 
of the apparent success of Model VI in the explanation 
of the X-ray properties of groups and clusters, it deserves 
further investigation of whether this disagreement arises 
from other effects. First, we exclude the contribution of 
nearby {z < 0.2) and massive (M > 5 x lO^^M©) clus- 
ters to the XRB, in the sense that these bright, diffuse 
X-ray sources may have been resolved by current deep X- 
ray observations. This leads to a moderate decrease of the 
expected XRB in high energy band but has only a minor 
effect on the soft XRB. Second, we attempt to include the 
absorption of X-ray emission by the neutral hydrogen in 
groups and clusters as a consequence of radiative cooling. 
To this end, we assiimc a simple King model for the den- 
sity distribution of the cooled material, in which we take 
the core size to be rc = O.lrvir- Moreover, we normalize 
this density profile using the total mass of the cooled IGM 
given by cooling scenario. Equation (14). The total X-ray 
luminosity at frequency u with self- absorption for a given 
halo reads 



dLx 
dhpV 



= 1.15 X 10' 



,56 



keVs- 
keV 



Jo y^ej 



Mpc/ ^0 



f^e J \ 10~2^ergs s~^ cm^ keV 

2 



(23) 



where a is the effective absorption cross-section, Ni is the 
hydrogen column density, and 



aNi = 0.31 



lO^^^cm- 



«H0 



10- 



r cos( 



r^i^ — sin^ 




in which nno is the central number density of neutral hy- 
drogen, r and rc are in units of Mpc. We adopt the fitting 
formula of Morrison & McCammon (1983) to calculate a. 
The XRB with the correction for hydrogen absorption is 
shown in Figure 5d. It appears that the inclusion of self- 
absorption results in a significant decrease of the XRB 
below E ~ 0.2 keV, while the rest XRB at £ > 0.2 keV 
remains almost imchanged. 

As for Model VII, the fact that the Lx-T relation at the 
low temperature slightly exceeds the observed data (see 
Fig. 3), despite the excellent agreement between the model- 
predicted entropy distribution for AS = 50-100 keV cm^ 
and the observations, indicates that the expected XRB 
at soft energy band will exhibit an excess relative to the 
observational limits. Indeed, although raising the extra 
entropy floor to AS = 100 keV cm^ leads to a decrease 
of the expected XRB relative to the cooling model predic- 
tion, in the low energy range E < 1 keV the prediction by 
model VII and observation limits still differ by a factor of 
~2. 

We notice, however, that the XRB spectra predicted by 
cooling, preheating and cooling plus heating models ac- 
tually have similar shape except that they are displaced 
upward in amplitude by a factor of 2-5 as compared with 
the unresolved soft XRB. This may suggest a common 
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origin of the discrepancies if they are not dne to our in- 
complete knowledge of gas physics in groups and clusters. 
We thus recalculate the XRB for all the models adopting 
a lower value of the normalization of the mass function of 
clusters and groups (jg = 0.7, instead of cts = 0.9, (e.g. 
Seljak 2001; Schuecker et al. 2002; etc.). The correspond- 
ing XRB spectra are illustrated in Figure 6. It appears 
that cooling model is still insufficient to bring the expected 
XRB to the observational limits in the soft energy band 
below E < 1 kcV. However, both preheating (especially 
Model V) and cooling plus heating models now yield the 
XRB spectra which agree nicely with the unresolved XRB. 
Perhaps, this indicates that the diffuse XRB could also be 
used as an independent constraint on the normalization of 
group and cluster abundance. 

5. CONCLUSIONS 

We have explored a set of simple analytic models for the 
distribution and evolution of the IGM in groups and clus- 
ters, aiming at highlighting the dominant physical process 
for the hot IGM, which may complement our understand- 
ing of the essential physics in addition to employment of 
hydrodynamical simulations. In particular, we have ad- 
dressed the question of whether the prevailing scenarios, 
namely, radiative cooling, preheating and a combination of 
cooling and heating, can explain both the global observed 
X-ray properties (entropy distribution and Lx-T relation) 
of groups and clusters and the observational limits on the 
contribution of the diffuse IGM in virialized halos to the 
XRB within the framework of standard CDM structure 
formation with an amplitude of matter power spectrum 
as = 0.9. Our main conclusions are summarized as fol- 
lows: 

Without radiative cooling and extra heating in addi- 
tion to gravitational shocks and adiabatic compression, the 
IGM appears to be too concentrated in halos to explain 
all the X-ray observations, especially in low-mass systems. 

Inclusion of radiative cooling may allow one to 
marginally reproduce the entropy floor seen in the cen- 
tral regions of groups and clusters (Voit & Bryan 2001b). 
However, it is still insufficient to account for the steepen- 
ing of the Lx-T relation on group scale, if X-ray surface- 
brightness bias is not included, and the observational up- 
per limits on the diffuse XRB. For the latter, the difference 
is by a factor of 2-5, consistent with the previous findings 
of Wu et al. (2001). Although cooling is certainly an im- 
portant process in the formation and evolution of galaxies, 
groups and clusters, and also for the explanation of the en- 
tropy floors shown in Fig.l, and even responsible for the 
scale-dcpendcncc of the IGM mass fraction (Wu & Xue 
2002b), energy feedback from galaxy formation should be 
included in order to resolve the inefficiency problem and 
the cooling crisis (Balogh et al. 2001). 

A phenomenological treatment of preheating, regardless 
of whatever the heating sources and mechanisms would 
be, may allow us to both reproduce the X-ray luminosity 
distribution of groups and clusters and explain the obser- 
vational limits on the XRB from IGM, provided that the 
level of entropy fioor can be raised to as large as 200 keV 
cm^ at present epoch. Apparently, the entropy profile in 
this case exceeds the measurements (Fig.lc). When we 
fix the value of input entropy floor to AS' = 120 keV cm^, 



which gives a nice fit to the observed central entropy distri- 
bution of groups and clusters, our predicted Lx-T relation 
and XRB both show an excess at low-temperature/energy 
range. In a word, a naive preheating model may meet dif- 
ficulty in the recovery of all the observational phenomena. 

Our cooling model with extra heating also fails when 
the predicted soft XRB is compared with the observa- 
tional limits. It has been shown that the inclusion of self- 
absorption by the cold gas in groups and clusters and the 
exclusion of the contribution of nearby, massive clusters 
to the XRB do not alter the conclusion significantly. Of 
course, our cooling plus heating scenario differs from the 
commonly used preheating plus cooling model in which the 
IGM was preheated before radiative cooling comes into ef- 
fect. 

Despite of their oversimplification, our analyses as a 
whole suggest that none of the present models can si- 
multaneously account for the observed X-ray properties 
of groups and clusters and the residual soft XRB. Actu- 
ally, the latter constitutes a very stringent constraint on 
the content and physical processes of the IGM in groups 
and clusters. Recall that the current observational limits 
even do not account for the fact that a large fraction of 
the unresolved soft XRB (~ 10% of the total XRB) which 
we have adopted in this paper may still arise from the X- 
ray emission of discrete sources (e.g. faint galaxies). More 
sophisticated models in combination with hydrodynami- 
cal simulations will thus be needed to further address the 
issue. 

Finally, the theoretically predicted XRB depends sen- 
sitively on the normalization of group and cluster abun- 
dance. If a lower value of (Tg = 0.7 is adopted instead of 
(7s = 0.9, it has been shown that the expected XRB spec- 
tra from preheating model and cooling plus heating model 
agree nicely with current observational limits. In other 
words, the observational constraints on the unresolved dif- 
fuse XRB lends support to a lower amplitude of matter 
power spectrum as « 0.7 (e.g. Seljak 2001; Schuecker 
et al. 2002). Indeed, the discrepancy between theoretical 
predictions (at least for preheating model and cooling plus 
heating model) and X-ray observations, if it is not a re- 
sult of the oversimplification of our analytic models, may 
have simply arisen from the too high normalization of the 
present matter power spectrum. 

We thank an anonymous referee for valuable sugges- 
tions. This work was supported by the National Sci- 
ence Foundation of China, and the Ministry of Science 
and Technology of China, under Grant No. NKBRSF 
G19990754. 
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Fig. 1. — Entropy distributions of groups and clusters measured at O.lrvir- Observational data are taken from Ponman et al. (1999; PCN), 
Lloyd-Davies et al. (2000; LPC) and Xu et al. (2001; XJW). Dash-dot-dot-dot line is the self-similar model (Model I) prediction. For cooling 
model (a) the entropy profile at 2 = 1 is also illustrated. 
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Fig. 2. — The X-ray luminosity - temperature relations predicted by cooling models II and III are compared with observations, in which 
we consider a time-evolving metallicity (upper panel) and a constant metallicity of 0.3Zq (lower panel). Meanwhile, in the upper panel we 
demonstrate how the Lx-T relation is modified by different X-ray flux limits Siimit in units of ergs s~^ arcmin"'^ cm~^ (grey lines). The 
Ly^-T relations at redshift 2 = 1 for the two models are also illustrated (dashed lines). For comparison, the results for self-similar model are 
shown by the dash-dot-dot-dot lines. 
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Fig. 3. — The X-ray luminosity - temperature relations predicted by preheating models IV (upper panel) and V (lower panel). For Model 
IV we show the results for two energy inputs: kAT = 0.25 and 0.3 keV, respectively. Note that there is a slight excess of the predicted Lx 
relative to the observed data in the temperature range from ~ 2 to ~ 8 keV. For Model V two different constant entropy floors are added to 
the IGM: AS = 120 and 200 keV cm^, in corresponding to the entropy distributions in Fig.l. For both Models IV and V, the Ly^-T relations 
at 2 = 1 are also shown (dashed lines). 
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Fig. 4. — The X-ray luminosity - temperature relations predicted by cooling plus heating models VI (upper panel) and VII (lower panel). 
Dashed lines are the Ly^-T relations at z = 1. 
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Fig. 5. — Comparison of the predicted XRB with the observational upper hmits on the contribution of diffuse IGM to XRB. Grey Unes: 
the self-similar results with time-varying (solid) and constant (dashed) metallicity. Downarrows are the observational upper limits in the 2-10 
keV band (Mushotzky et al. 2000; Giacconi et al. 2001; Tozzi et al. 2001). Filled triangles are the constraints obtained by Cowie et al. (2002) 
in the 2-8 keV band. Crosses, open squares and open circles in the 1-2 keV range correspond to three evaluations of the total XRB intensity 
by Chen et al. (1997), Miyaji et al. (1998) and Gendreau (1997), respectively, after the source contributions detected by Hasinger et al. 
(1998), Mushotzky et al. (2000), Giacconi et al. (2001) & Tozzi et al. (2001) are removed. Filled circles and diamonds are the upper limits 
derived by Bryan & Voit (2001) in the 0.5-2 and 0.1-0.4 keV bands, respectively. The total XRB estimated by summing up the contributions 
of groups and clusters characterized by their X-ray luminosity functions and the observationally determined Ly^-T relations is shown by the 
shaded region (Wu & Xue 2001). 
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Fig. 6. — The same as Fig. 5 but for as = 0.7. 



